Relaxation of a Colloidal Particle into a Nonequilibrium Steady State 
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We study the relaxation of a single colloidal sphere which is periodically driven between two 
nonequilibrium steady states. Experimentally, this is achieved by driving the particle along a toroidal 
trap imposed by scanned optical tweezers. We find that the relaxation time after which the probabil- 
ity distributions have been relaxed is identical to the decay of the velocity autocorrelation function, 
measured in a steady state. In quantitative agreement with theoretical calculations the relaxation 
time strongly increases when driving the system further away from thermal equilibrium. 
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The understanding of thermodynamic processes at 
small length scales is of central importance at the inter- 
face of physics, biology and chemistry. Classical thermo- 
dynamics, as originally developed for macroscopic sys- 
tems with many internal degrees of freedom cannot be 
applied to e.g. molecular machines, proteins or micro- 
mechanical devices. This is because at microscopic scales, 
thermal fluctuations must not be neglected and the fa- 
miliar well-defined thermodynamical quantities have to 
be replaced by corresponding distributions of finite width 
[H, 0,13 ■ The situation is further complicated when these 
systems are driven out of thermal equilibrium as often en- 
countered within their natural environment. The treat- 
ment of fluctuations in such nonequilibrium situations is 
even more difficult since it requires the full knowledge 
of the system's dynamics. Despite considerable progress 
in deriving exact relationships which are valid beyond 
thermal equilibrium [1, 0, Q a comprehensive theoretical 
description of nonequilibrium is still lacking. 

Among the huge manifold of nonequilibrium condi- 
tions, nonequilibrium steady states (NESS) are certainly 
the most simple conceivable situations, being character- 
ized by a time-independent probability distribution in 
the presence of a non- vanishing probability current. Ac- 
cordingly, NESS present ideal conditions for fundamen- 
tal studies and tests of nonequilibrium properties on 
a microscopic scale 0, . 

In this Letter we experimentally investigate the re- 
laxation behavior of a single colloidal particle which is 
periodically driven between two different nonequilibrium 
steady states NESS I and NESS II being created by scan- 
ning optical tweezers. We find that the NESS relaxation 
time as defined by the decay of the probability distri- 
bution only depends on the final state but is indepen- 
dent of the initial one. In addition, we show that this 
relaxation time is identical to that obtained by the de- 
cay of the velocity autocorrelation function in the steady 
state regime, i.e. after relaxation has been completed. 
In agreement with theoretical calculations, the relaxation 
time increases when driving the system further away from 
thermal equilibrium. 

The experimental setup has been already described 
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FIG. 1: (color online) (a) Realization principle of the cre- 
ation of a NESS for a colloidal particle by scanning a focussed 
laser beam, (b) Schematic representation of NESS I and 
NESS II which correspond to a tilted periodic potential. A 
sudden change of the driving force / and phase (f> = n leads to 
a redistribution of the related probability distributions pi and 
Pii • 



elsewhere and will be discussed here only in brief [lo| . 
Well-defined nonequilibrium steady states for a colloidal 
silica particle immersed in water, with radius a = 
0.65 //m are created by scanning the highly focussed 
beam of a Nd:YAG laser (A = 532 nm) along a circle with 
radius R = 1.14 /xm (see Fig. [1]). At rather high scan- 
ning frequencies the particle cannot follow the tweezers 
motion due to the viscous forces of the fluid, and it is 
confined to an effective three-dimensional toroidal opti- 
cal trap. At intermediate scanning frequencies, however, 
each time the scanning laser focus passes the particle, 
a small displacement of the colloid along the scanning 
direction is induced. Since individual kicks are not re- 
solvable by digital video microscopy 



2l|, in this regime 



the scanning tweezers can be considered as exerting a 
constant force / on the particle along the angular co- 
ordinate X 



12| . For a scanning frequency of 
200 Hz and a laser intensity /o ~ 40 mW this leads to 
a drift velocity oi v k, 7 /im/s. In addition, the laser in- 
tensity is weakly modulated along the toroidal trap with 
an electro-optical device whose input signal is synchro- 
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FIG. 2: (color online) (a) Particle trajectory (solid line) and 
p{x,t) as gray scaled background. Dashed vertical lines in- 
dicate the transitions between two different NESS with du- 
rations ti = 5 s, tii = 2 s and a phase difference 4)u — 
TT. (b) Normalized steady state probability distributions of 
NESS I (open bars) and NESS II (closed bars). 
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FIG. 3: (color online) (a),(c) Open symbols: Mean drift ve- 
locity {v{t)) after switching from NESS II to NESS I and vice 
versa (/i 53 keT/^im, 0i = 0). The NESS II parameters 
are given m Tab. |T](a): (i), (c): (iv). Solid line; Exponen- 
tially damped sinusoidal function. (b),(d): Corresponding 
potentials (inset) and steady state distributions of NESS II . 



nized with the scanning motion of the laser focus. For a 
periodic intensity modulation I(x) = /o + A/sin(x) this 
leads to additional optical gradient forces, i.e. a static 
potential V{x) = —-^ am{x + (j)) acting on the particle. 
The value of can be controlled by the relative phase 
difference between the scanned tweezers motion and its 
intensity variation. In total, the colloid is subjected to a 
tilted periodic potential U {x) = V{x) — fRx correspond- 
ing to a NESS where / and V{x) can be tuned by /g 
and AJ, respectively. The driving force / and V{x) are 
not known a priori but can be reconstructed via a gen- 
eralized Boltzmann factor from the measured stationary 
probability distribution p"^ (x) and the probability current 
in the system [ist - 

The relaxation of a colloidal particle into a NESS is in- 
vestigated by periodically toggling between two differing 
steady states. This is accomplished by a sudden change 
in the driving force / and the phase </> according to the 
protocol 

if < t < ti : /i , 01 = (NESS I ) 

if<i <i<<i +tu: (NESS II). ^ ' 

Unless otherwise stated, Vq is kept constant at Vq ~ 
100 keT. The duration times ti and are chosen suffi- 
ciently long to allow the system to reach the correspond- 
ing stationary probability distributions pf and pjj . The 
entire protocol is typically repeated up to 800 times dur- 
ing each experiment to obtain adequate statistical aver- 
ages. 

To illustrate the principle of our experiments, we first 
discuss the situation where NESS II is close to thermal 
equilibrium. This is achieved by applying a rather weak 
driving force fu ~ 4 keT/zim. Accordingly, U(x) ex- 
hibits a potential well of about 80 keT where the particle 
remains strongly localized, thus closely resembling equi- 
librium conditions (locked state). This is clearly seen by 
the trajectory (solid line) in Fig. [D^a) which is confined 



to a small range of x-values. In contrast, NESS I has a 
much stronger force fi « 53 keT/zim. Therefore, U{x) 
exhibits no local minimum and the particle is free to drift 
along the entire torus (running state). 

Because of the superimposed Brownian motion, the 
particle trajectory varies between each cycle of the proto- 
col. These fluctuations are taken into account by consid- 
ering the probability distribution p(x,t), i.e. the prob- 
ability of finding the particle at time t at position x. 
The measured p{x, t) obtained from about 800 cycles 
of the protocol is shown as gray scaled background of 
Fig. [2ja). For — 1 < t < s the particle has relaxed 
to NESS II where it is localized inside the deep poten- 
tial minimum. The corresponding strongly peaked steady 
state probability distribution pfj is shown as closed bars 
in Fig. [DJb). Upon suddenly switching to NESS I at 
t = s, the particle starts to circulate along the en- 
tire toroidal trap; this leads to a broadening and a shift 
of the maximum in p{x,t). The damped oscillatory be- 
havior of p{x ~ const., t) is typical for the relaxation into 
a nonequilibrium steady state. This is in contrast to the 
situation at t — 5 s when the protocol switches back to 
the equilibrium-like conditions of NESS II . Here, p{x = 
const., t) monotonically approaches its final value pfj . It 
should be noted that the relaxation from NESS I into 
NESS II proceeds much more rapidly than into the other 
direction. 

To quantify our findings we calculate the mean drift 
velocity {v{t)) = ^ ^(*+'^*)~^(*~^*) )^ which is obtained by 
averaging the actual particle velocity v{t) over several 
hundred cycles of the protocol. Since the length scale 
over which the potential U{x) varies is more than one 
order of magnitude larger than the maximal particle dis- 
placement between two consecutive {At = 33 ms) video 
frames, {v{t)) can be obtained from the experimentally 
determined trajectories. The symbols in Fig. EJa) show 
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1.4 ± 0.2 


(ii) 


4 








1.6 ±0.25 


(iii) 


38 


TT 


3.4 


1.6 ±0.25 


(iv) 


52 


TT 


9.7 


1.9 ±0.3 


(v) 


99 


TT 


21 


1.6 ±0.25 



TABLE I: Measured relaxation time of NESS I for different 
parameters of the initial NESS II . 




{v{t)) for the same data set as in Fig. [21 After switch- 
ing to NESS I , {v{t)) is a decaying oscihatory function 
which converges to the corresponding mean steady state 
velocity. For < > 1 s it can be weh described by an ex- 
ponentially damped sinusoidal function (solid line) with 
decay time r = 1.4±0.2 s and the oscillation period given 
by the mean particle revolution time Tn = 0.9 s. Similar 
as above, {v{t)) for the relaxation into the equilibrium- 
like NESS II is purely exponential with a decay time of 
0.3 s. 

In order to understand how the relaxation into a NESS 
compares with that into thermal equilibrium we consider 
the relaxation time t°'^ of an overdamped Brownian par- 
ticle into a parabolic potential. According to Ornstein- 
Uhlenbeck [IJ, [l5 1 this is given by 

Qirrja 



(2) 



where rj is the viscosity of the solvent, k the potential 
curvature and a the particle radius. Obviously, t^"^ is en- 
tirely determined by the state into which the relaxation 
occurs and independent of the initial conditions. In order 
to investigate whether this holds also for relaxation pro- 
cesses into noncquilibrium states, we systematically vary 
the initial NESS II (by changing /n and 0//) and study 
the relaxation into the identical final NESS I (same pa- 
rameters as in Fig.[2|). As an example, Fig.[3^c) shows the 
relaxation for /n ~ 52 keT/zim and 0// = tt. It should 
be realized that although the steady state distribution of 
the initial state in Figs. [3Jb) and (d) is rather different, 
the decay of the mean drift velocity is - within our exper- 
imental errors being caused by the finite number of tra- 
jectories and small optical drifts - identical. This is also 
seen m Table U which summarizes five relaxation experi- 
ments from different NESS II into the identical NESS I . 
Within our experimental accuracy we observe the same 
relaxation time = 1.6 ± 0.2 s. At least in case of 
the specific NESS as considered here, this suggests that 
the relaxation time only depends on the final state. We 
confirmed the independence of the relaxation time for a 
variety of different NESS I conditions. Due to techni- 
cal details all these experiments were performed in the 
running regime. 

According to the fluctuation dissipation theorem 
(FDT), the temporal decay of fluctuations does not de- 
pend on whether they arc imposed by an external force 



FIG. 4: (color online) (a) Measured velocity autocorrelation 
function (symbols). The solid line is an exponentially decay- 
ing sinusoidal fit. (b) Relaxation time vs. Tcorr for different 
NESS where the driving force has been varied between 40 and 
87 ksT/fJ^ra. The open and closed symbols correspond to po- 
tential depths of 125 and 100 keT, respectively. The straight 
line has slope one. 



or spontaneously generated by the system itself. How- 
ever, it is important to realize that this identity is only 
valid in or close to thermal equilibrium [l6| . There- 
fore it is not a priori clear whether in a driven sys- 
tem (as considered here) , i.e. the decay time in re- 
sponse to a sudden change of the external driving force 
is identical with the decay of the steady state fluctua- 
tions. As shown in Fig. IH^a) the mean drift velocity au- 
tocorrelation function C{t) = {v{t')v{t' + t))t', obtained 
via a stationary measurement under NESS I conditions, 
is an exponentially decaying sinusoidal function with a 
decay time of Tcon- = 1.7 ± 0.2 s. Within the ex- 
perimental error this value is again identical with the 
above determined . To test whether this agreement is 
generally valid, we performed additional measurements 
with different driving forces / and potential depths Vq. 
In Fig. [3]Jb) we compare Tcorr and Tr as measured for 
40 kBT//im < // < 87 keT/^m and Vo = 100 keT 
and Vo = 125 keT, respectively. The good agreement 
between the data points and the solid line (slope one) 
supports that Tr = Tcorr ^nd suggcsts that (as in equi- 
librium), the relaxation time of a NESS can be measured 
via transient or stationary measurements. 

In order to compare the relaxation time with theory, 
we calculate Tr by numerically solving the Fokker-Planck 
equation 



dtp{x,t) = -dx[ldoF{x) - Doda:]p{x,t), 



(3) 



with R the torus radius and F{x) = the total 

external force acting on the particle. The transport coef- 
ficients are assumed to be uneffected by the external driv- 
ing force [1, 0] , therefore the free diffusion coefficient Dq 
and the mobility fiQ are taken from thermal equilibrium. 
In units of dimensionless time i = {Do/R'^)t = e~^t and 
force F{x) = {R/kBT)F{x) the Fokker-Planck equation 
reduces to dtp{x, i) — Lxp{x, i). Since the Fokker-Planck 
operator = —dxF{x) + 9^ has no explicit time depen- 
dence, a separation ansatz for the probability distribu- 
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FIG. 5: (color online) Symbols: Measured decay times 
Tcorr for two different driving forces. The solid lines show 
the parameter- free numerical prediction for eri. For poten- 
tial depths smaller than the critical amplitude Vc ~ 'i.fR the 
minimum in the tilted potential vanishes, which leads to an 
enhancement of the diffusion coefficient 

iion p{x,i) ~ exp(— A„ t)g„(a;) leads to the following 
eigenvalue equation 

- A„ qn {x) = q„ (x) . (4) 

The relaxation of an arbitrary given initial probability 
distribution is described by the complete set of eigenval- 
ues A„. However, in the long time limit only the two 
smallest eigenvalues Aq and Ai are relevant. The station- 
ary solution p^{x) is given by Aq = and go (2;). The 
real part of Ai, TZ{\i) = t^^ determines the asymptotic 
time dependence of the relaxation process. Therefore the 
relaxation time is £Ti . Since Eq. ^ has no analytical so- 
lution, for the determination of the eigenvalues we have 
to expand the eigenfunctions into an orthonormal ba- 
sis. Due to the periodic nature of the system a suitable 
choice is the Fourier series qn{x) = c|"'' exp(i^x). 

A straightforward calculation leads to 

-A„4")=5]L,.,cJ"\ (5) 
I 

an eigenvalue equation for the matrix L = (Lj.; ) . In case 
of the experimentally realized sinusoidal potential, L is 
tridiagonal [l3| • After truncating the size of the matrix 
to a finite value its eigenvalues are easily found using 
standard numerical algorithms. 

Fig. [5] shows the calculated en (sohd lines) as a func- 
tion of the potential depth Vb and for two different driv- 
ing forces /. For large Vb the system becomes equilibrium 
like and the relaxation time asymptotically approaches 
the Ornstein-Uhlenbeck result of Eq. ([2]). At small Vb 
the relaxation time is only determined by the timescale 
the particle needs to diffuse along the toroidal trap, i.e. 
STi = Dq/R^. These two limiting cases are connected 
via a monotonic curve. The closed symbols correspond 
to the experimentally determined decay time Tcorr • The 



excellent parameter free agreement between experimen- 
tal and numerical results again supports our assumption 
that Tcorr IS cqual to Tr (see Fig. |4]) and a posteriori jus- 
tifies that even for driven colloidal systems the Fokker- 
Planck equation, with Dq taken from equilibrium, is still 
valid. 

In summary, we have investigated the relaxation be- 
havior of a colloidal particle into a NESS. Our results 
show that the NESS relaxation time is independent of 
the initial conditions from which the relaxation process 
starts. In agreement with calculations we confirm that, 
in case of a driven colloidal particle, the noncquilibrium 
relaxation time is identical to the decay time of the veloc- 
ity autocorrelation function. It must be emphasized that 
it is not clear yet whether our observations are generally 
valid to arbitrary NESS or restricted to particular situa- 
tions. We hope that our work will stimulate further theo- 
retical studies in this direction. It will be also interesting 
to perform similar stationary noncquilibrium relaxation 
measurements in systems of sheared polymers [lij or vesi- 
cles 0. 
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